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Abstract 


A  Fourier-integral  method  is  developed  to  obtain  transient  solutions  to  potential  wavemaker 
problems.  This  method  yields  solutions  that  are  uniformly  valid  for  wavemaker  velocities  which 
need  not  be  given  as  powers  of  time.  The  results  are  compared  with  known  small-time  and  local 
solutions.  Examples  considered  include  ramp,  step  and  harmonic  wavemaker  velocities.  As  time 
becomes  large,  the  behavior  near  the  wave  front  is  derived  for  the  impulsive  wavemaker,  and  for 
the  harmonic  wavemaker  it  is  shown  that  the  steady-state  solution  is  recovered.  The  soluticm  for  a 
wavemaker  velocity  given  as  a  Fourier  cosine  series  compares  favorably  with  the  computational 
and  experimental  results  of  Dommermuth  et  al.  (1988).  Capillary  effects  are  included  and 
nonlinear  effects  are  discussed. 
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1  Introduction 


The  study  of  water  waves  associated  with  surface-piercing  bodies  has  long  been  an  interesting 
and  important  area  in  fluid  mechanics.  In  many  cases,  however,  both  analysis  and  computation 
run  into  difficulties  near  the  intersection  of  a  body  with  the  free  surface.  To  examine  this  problem, 
a  number  of  researchers  have  considered  a  vertical  wavemaker  moving  horisontally  in  a  fluid  of 
finite  depth.  Peregrine  (1972)  used  a  moving  coordinate  system  attached  to  the  wavemaker  in  his 
small-time  expansion,  noted  a  singular  behavior  at  the  contact  line,  and  explained  the  necessity 
of  a  local  solution.  Chwang  (1983)  obtained  a  s<Jutioa  with  a  stationary  coordinate  system  and 
argued  that  the  singularity  lies  outside  of  the  physical  domain.  However,  he  did  not  properly 
expand  the  moving  wavemaker  boundary  conditions  about  the  singular  point.  The  logarithmic 
singularity  in  the  small-time  solution  was  confirmed  by  Lin  (1984)  by  a  Lagrangian  description 
of  the  problem. 

It  was  not  until  the  work  of  Roberta  (1987)  that  a  successful  local  solution  for  the  power- 
law  movement  of  the  wavemaker  was  obtained  for  small  time  and  small  Firoude  number.  He 
found  that  for  small  time  the  solution  varies  significantly  in  the  neighborhood  the  contact  line 
and  gave  a  self-similar  formulation  to  describe  this  behavior.  If  the  wavemaker  starts  to  move 
impulsively  (“step  velocity"),  however,  the  ne^ected  nonlinear  effect  becomes  important  close 
to  the  wavemaker,  since  the  linearised  vertical  velocity  becomes  infinitely  large  as  time  becomes 
smaller. 

The  present  study  ^proaches  the  problem  using  inviscid-flow  theory  and  develops  a  simple 
method  to  obtain  solutions  that  are  uniformly  valid  for  wavemaker  velocities  which  need  not  be 
given  as  powers  of  time.  It  is  also  possible  to  include  the  effects  of  surface  tension  and  initial 
free-surface  elevation  caused  by  non-sero  static  contact  an|^e  with  little  added  effort.  For  viscous 
fluids  the  unsteady  motkm  of  this  contact  line,  with  the  no-slip  condition  incorporated  at  the 
body,  would  requite  special  consideration,  as  discussed  in  the  review  by  Dussan  V.  (1979)  and  in 
an  application  to  water  waves  by  Hocking  (1987).  We  will  not  consider  such  complice* here. 
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The  singularity  at  the  contact  line  can  be  shown  to  be  a  consequence  of  an  improper  expansion. 
The  small-time  s<^ution  ceases  to  be  a  valid  asymptotic  series  in  time  as  the  distance  from  the 
contact  line  decreases  because  the  higher-order  terms  do  not  remain  small.  A  correct  local 
solution  requires  a  new  length  scale  that  varies  in  time.  In  other  words,  very  near  the  contact 
line  spatial  variables  are  coupled  with  time,  so  that  the  small-time  expansion  does  not  yield  a 
well-behaved  solution  there.  An  alternative  to  the  small-time  expansion  is  an  expansion  for  small 
Froude  number  followed  by  a  Laplace  transformation,  as  used  by  Roberts  (1987),  at  the  Fourier 
transformation  used  in  the  present  work.  Both  of  these  methods  can  be  applied  to  obtain  a 
solution,  since  neither  requires  separation  of  variables.  However,  the  latter  can  be  more  powerful 
and  versatile  in  that  the  interpretation  of  the  results  and  the  extension  to  arbitrary  wavemaker 
velocities  are  accomplished  more  easily. 

We  begin  in  §2  by  introducing  the  fommlation  of  the  vertical  wavemaker  problem.  The 
velocity  of  the  wavemaker  is  given  as  a  function  of  time,  and  the  surface  tension  on  the  free 
surface  is  retained  with  non-aero  static  contact  angle. 

In  §3  and  §4  we  consider  two  hypothetical  wavemaker  velocities,  expressed  by  ramp  and 
step  functions  in  time.  The  solutions  for  sero  surface  tension  are  shown  to  agree  with  the  local 
solutions  of  Roberts  (1987)  for  small  time  near  the  contact  line.  Also,  sufficiently  far  from  the 
contact  line,  they  are  shown  to  agree  with  the  small-time  solution  of  Peregrine  (1972).  The 
nonlinear  fmnulation  and  large-time  behavior  for  the  step  velocity  are  also  discussed. 

More  general  types  of  wavemaker  velocities  are  discussed  in  $5  and  §6.  One  example  is  a 
noore  realistic  velocity  that  starts  from  lero  and  increases  toward  a  finite  constant  value.  It  is 
shown  how  the  solutions  for  step  and  ramp  velocities  can  be  recovered  as  limiting  cases.  In  §6, 
the  present  method  allows  a  transient  solution  for  a  simple-harmonic  wavemaker.  Sufficiently  far 
from  the  contact  line  and  in  the  limit  as  time  approaches  infinity,  the  soluti<m  agrees  with  that 
of  Havelock  (1929),  as  cited  by  Yih  (1979).  As  a  final  example  in  §6,  we  examine  a  wavemaker 
velocity  considered  by  Dommennuth  et  al.  (1988)  and  compare  our  analytical  solution  with  their 
computatkmtU  and  experimental  results. 
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This  study  is  prompted,  in  part,  by  computational  difficulties  caused  by  bodies  intersect¬ 
ing  free  surfaces.  These  include  small-wavenumber  oscillations  close  to  the  contact  line  that 
may  imply  a  physical  or  numerical  instability.  Normally,  the  spatial  dependence  is  solved  by  a 
boundary-integral  approach,  but  these  techniques  are  known  to  have  difficulties  with  corners, 
even  when  singular  behavior  is  not  present  (Schults  tc  Hong,  1988).  The  standard  approach  is  to 
separate  the  spatial  and  temporal  behavior,  thus  introducing  singular  behavior  into  the  problem 
as  mentioned  above,  and  it  is  only  the  inability  of  the  boundary-integral  computations  to  model 
this  apparent  singularity  that  allows  the  excellent  agreement  (at  least  for  regions  not  too  close 
to  the  wavemaker)  of  previous  computations  (Dommermuth  et  al.,  1988). 

2  Formulation 


We  consider  the  fluid  motion  due  to  a,  moving  wall  as  shown  in  Fig.  1.  If  the  fluid  is  inviscid 
and  incompressible,  and  the  motion  starts  from  rest,  the  flow  will  be  irrotational  according  to 
Kelvin’s  theorem  and  is  described  by  the  Laplace  equation.  In  nondimensional  variables 


=  0  for  ar  >  t(t),  -1  <  y  <  f/(»,  <),  (1) 

where  ^(z,y,t)  is  the  velocity  potential,  *;(s,t)  is  the  free-surface  elevation  measured  from  the 
undisturbed  level  at  infinity,  and  s(t)  is  the  dispiac«nent  of  the  wall  from  its  initial  location, 
which,  of  course,  is  the  time  integral  of  the  wavemaker  velocity.  The  velocity,  length,  and  time 
scales  are  chosen  to  be  h,  y/h/g,  where  h  is  the  undisturbed  depth  of  the  fluid  and  g  is  the 
gravitational  acceleration.  The  fluid  velocity  at  the  wavemaker  is  prescribed  as 


=  o;«(t)  on  *  =  s(t), 

where  a  is  the  Fkoude  number.  For  example,  if  the  dimensional  velocity,  17,  is  given  by 

<>0, 


where  C  and  q  are  constants,  we  have 


a  = 


(2) 

(3) 
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Thus,  u(<)  and  ${t)  in  (2)  become 


«(<)  =  «•,  s(0  =  at«+V(«+l). 

In  paiticular,  the  cases  for  9  =  0  and  9  =  1  correspond  to  a  step  velocity  and  a  ramp  velocity, 
respectively. 

The  aforementioned  u(t)  describes  a  power-law  movement  of  the  wavemaket,  which  can  ap¬ 
proximate  the  earlier  stages  of  a  more  general  motion  for  small  time  and,  thus,  is  frequently 
used  in  small-time  analyses.  In  the  present  work,  however,  the  general  expression  for  u(t)  will 
be  maintained  and  a  solution  valid  for  all  time  will  be  sought. 

In  the  presence  of  non-sero  surface  tension,  the  kinematic  and  dynamic  boundary  conditions 
on  the  free  surface  become 


^  on  y  =  1) 

where  the  nondimensional  surface  tension  T,  the  reciprocal  of  the  Bond  number,  is  defined 


T  = 


pgh*' 


(4) 

(5) 

by 


Here,  a  is  the  surface-tension  constant,  and  p  is  the  density  of  the  fluid.  On  the  bottom,  the 
vertical  velocity  component  vanishes,  so  that 


dr  ~  0  OD  y  —  — 1- 

Since  the  motion  starts  from  rest,  the  initial  conditions  are 


(6) 


1  <0 

(7) 

t  <  0. 

(8) 

Here,  a  is  a  constant  determined  by  a  =  tan(r/2  -  9,),  where  9,  is  the  static  contact  angle. 
The  initial  free-sutfare  ekvstion  (8)  satisfies  the  static  linear  equivalent  of  (5)  for  small  a  and 


becomes  sero  in  the  limit  as  the  surface  tensi(»  becomes  sero  or  the  initial  contact  angle  becomes 
90*. 

Instead  of  introducing  a  small-time  expansion,  which  is  not  unif(»mly  valid  near  the  contact 
line,  as  explained  in  $1,  we  now  use  an  expansion  for  small  Ftoude  number.  The  velocity  potential 
and  the  free-surface  elevation  are  expanded  as 


(*.  y.<)  +  y.  0  +  •  •  • 

»c>/5rexp(-^)  +  ariiix,  t)  +  o*ih(x.  <)  +  ••■. 


(9) 

(10) 


Although  T  is  typically  very  small,  it  is  possible  to  retain  T  =  0(1)  in  much  of  the  following 
with  very  little  added  complexity. 

Expanding  the  free-surface  boundary  conditions  and  the  boundary  condition  on  the  wave- 
maker  about  y  =  0  and  z  =  0,  respectively,  pves  the  equations  to  leading  <»der,  0(a),  as 


+  dipt  ~  9 

for  z  >  0,  -1  <  y  <  0 

(11) 

dig  =  u(t) 

on  z  =  0 

(12) 

dip  = 

on  y  s  0 

(13) 

Tmgg  =  0 

on  y  =  0 

(14) 

-a. 

aw 

II 

o 

1 

II 

§ 

(15) 

The  boundary  condition  (12)  requires  the  distance  s(f)  to  be  small,  so  that  a  restriction  should 
be  imposed  on  t,  depending  on  the  velocity  u(t).  For  the  step  velocity,  for  example,  the  condition 
is  t  (l/o)-  This  constraint  could  be  relaxed  by  applying  a  simple  coordinate  transformation 
z'  =  z— s(t)  to  fix  the  location  of  the  wavemaker  at  z'  =  0,  which  leaves  the  first-order  equations 
(11)-(15)  unchanged.  However,  for  large  time  it  should  be  expected  that  the  cumulative  effect  of 
omitted  nonlinear  terms  will  no  longer  be  negligible,  and  that  qj  and  ^  will  become  large.  F<» 
example,  siace  the  nmdimensional  wave  speed  in  shallow  water  is  dx/dt  :=  1  -b  0(a),  we  may 
antfcipstc  that  the  e]q>ansions  (9)  and  (10)  are  valid  for  large  t  only  if  at  is  small. 
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In  solving  (11)-(15),  we  first  decompose  two  psrts: 

«  1 

di  =  2ii(i)  2^  — e  ‘•'sininy  +  dKr.y.O.  (16) 

fixO  ** 

where  =  (n  -f  ^)t.  The  series  on  the  right-hsnd  side  satisfies  the  Laplace  equation  and  ail  the 
boundary  conditions  except  on  the  free  surface,  where  it  becomes  sero.  The  remaining  term  dt , 
then,  can  be  considered  as  a  correction  that  enables  the  cmnplete  solution  to  satisfy  appropriate 
free-surface  boundary  conditions.  Substitution  of  (16)  into  (11)-(15)  yields 


=  W 

foe  *  >  0,  -1  <  y  <  0 

(17) 

0 

II 

o 

II 

H 

§ 

(18) 

-|u(t)  ln(tanh  ^)  +  dj,  =  ifti 

on  y  =  0 

(19) 

~  Ttfigg  =  0 

on  y  =  0 

(20) 

0 

II 

•e- 

on  y=s— 1. 

(21) 

The  solution  for  is  sought  as  a  Fourier  cosine  integral: 


v4(fc,t)  cosh  *(y+ 1)  cos  Irxdi, 


which  already  satisfies  (17),  (18),  and  (21).  The  solution  for  171  is  then  also  a  Fourier  cosine 


integral: 


fh  =  f  B(k,t)  COB  kxdk. 
Jo 


Substituting  the  representations  (22)  and  (23)  into  the  free-surface  conditions  (19)  and  (20)  and 
eliminating  B  gives 


Att  cosh  k  +  k{l  +  Tk^)A  sinh  *  =  -^^(  J  +  Tk)  tanh  k,  (24) 

by  making  use  of 

y*"  tanh  fe  cos  = -b(taiih^)  (25) 

(Gradshteyn  k  Ryshik,  19M)).  Equation  (24)  has  a  general  solution 


A  =  i4^(k,t)-|-ct(k)sin^4-ea(k)cas/}(. 


(26) 
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where 


0=  y/k{l+Tk^)U^. 

Here,  i«  «  particular  solution  o(  (24)  for  a  given  u(t),  and  cj  and  ej  are  to  be  determined  from 
the  initial  conditions. 

The  Laplace  equation  does  not  requite  any  initial  condition.  Since  the  time  derivatives  appear 
only  in  the  free-surface  boundary  conditions,  the  initial  conditions  are  applied  on  the  free  surface. 
The  initial  condition  for  q,  (8),  is  converted  into  a  condition  for  d  through  (5).  Using  (9)-(10) 
again,  we  then  obtain 

s  0  and  dll  =  0 

at  y  =:  0  and  (  =  0.  Since  the  series  in  (16)  is  always  sero  at  y  =  0,  through  (22) 

yl(h.0)  =  ^(i.0)  =  0.  (27) 

The  solution  (26),  then,  is  completely  determined,  and  so  is  dl< 

Once  dt  >•  obtained,  the  first-order  free-surfsce  elevation,  iji,  is  given  either  by  (19)  or  by 
(20).  The  vertical  velocity, 

diy  =  -•-«(<)  ln(**®l>  0  •“!>  k  CKMkxdk, 

»  *  Jo 

has  a  logarithmic  singularity  in  the  first  term,  which  is  canceUed  by  the  same  singularity  in 
the  second  term.  This  will  be  examined  in  greater  detail  for  each  specific  case  in  the  fdlowing 
sections. 

The  highetHwder  velocity  potentials  also  satkiy  the  Laplace  equation  with  a  Neumann  c(»di- 
tion  <»  the  wavemaker,  so  that  the  same  Fburier-integral  method  as  for  di  can  be  used.  However, 
in  most  cases  we  do  not  evaluate  the  Fourier  integral  for  di  or  >71  to  obtain  exact  dosed-form 
solutions.  Thus,  the  nonhmnogeneous  terms  in  the  hi|^er-order  analysis  contain  smne  products 
of  Fourier  integrab,  which  makes  numerical  analysis  inevitable  except  for  certain  cases. 
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3  Ramp  velocity 


The  poirer-law  behavior  crfthe  wavemaker,  aa  given  by  (3),  deaervea  apecial  attention,  becauae 
it  can  expoae  the  initial  evolution  of  the  fluid  motion  without  unnecesaary  complication.  Due 
to  the  obvioua  diatinction  between  the  ramp  and  atep  velocitiea,  aeparate  conaideration  will  be 


given  to  each  caae. 

The  ramp  velocity  repreaenta  a  wavemaker  that  atarta  from  reat  and  increaaea  in  apeed  linearly 
with  time.  It  correapooda  to  the  caae  when  4  s  1  in  (3),  ao  that  the  Froude  number,  a,  and  the 
dimenaionleaa  velocity,  u(t),  in  (2)  beconae 


a  = 


C 

9 


and  w(l)  =  (. 


(28) 


Therefore,  the  expanaiona  in  (9)-(10)  require  amall  acceleration  of  the  wavemaker. 

Application  d  the  Fourier-integral  method  explained  in  the  previoua  aection  yielda  (24)  with 
'a(i)  =:  t.  The  aolution  aatiafying  (27)  ia 


A(k,t) 


2  1  / 
wk^eo^k  V  P  )' 


(29) 


Then,  from  (22)  and  (10),  the  flrat-order  velocity  potential  ia  obtained  aa 

coahk(y-t- 1) 


coahk 


coa*x||.  (30) 


The  vertical  velocity  of  the  fluid  on  the  free  aurface,  or  the  left-hand-aide  of  (19),  becomea 

2  ain^  ^  ,  dk  . 

di,  =  -y^  — ~  tanhkcoa**  — ,  (31) 

where  the  logarithmic  term  ariaing  from  the  infinite  aeriea  haa  been  cancelled  with  the  help  of 
(25).  The  free-aurface  elevation  obtained  by  uaing  (19)  ia 


when  «c  =  0. 

To  recover  the  amall-tiine  aolution,  we  can  rewrite  (32)  aa  the  aum  of  an  integral  from  sero 
to  K  and  an  integral  from  K  to  infinity,  with  K  choaen  auch  that  1  <  <  (!/(’)■  For 
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0  <  k  <  K,  CM0t  can  be  expanded  in  a  Tkylor  seriea  in  time,  whereas  tof  K  <  k  <  oo  we  have 
tanh  ib  =  1  +  0(e~*).  Then,  in  the  absence  of  surface  tensicm  (T  =  0),  (32)  can  be  expanded  as 
t  — >  0  with  z  fixed  to  give 

The  first  term  in  (33)  is  consistent  with  the  small-time  soluti(»  of  Peregrine  (1972)  and  Chwang 
(1983),  and  the  additional  terms  can  be  obtained  by  extending  the  small-time  expansicm  to 


o(an 

If  both  X  and  t  are  small,  with  x  =  O(t’),  the  integral  (32)  can  be  simplified  when  T  =  0  by 
the  addition  and  subtraction  of  tanh  k  in  the  numerator  of  the  integrand.  The  logarithmic 
term  is  obtained  explicitly  by  the  use  of  (25),  and  the  remaining  integral  is  written  in  two  parts 
as  before.  It  is  found  that 

—  cos  y/kt ' 


crt>  rx  2o<’  ^  A  l-cosVkA  dk 

X  ^  A  ~  r  Jo  \2~  it*  ^  k 


(34) 


(Joo  et  al.,  1988),  where  the  integral  is  a  function  of  c/t*.  Using  integration  by  parts,  we  find  the 
largest  term  to  be  f'*/(720x*)  as  s/(*  -*  oo,  which  agrees  with  the  expansion  of  (33)  as  x  -»  0. 
The  velocity  components  can  be  found  in  a  similar  way.  In  complex  form,  with  x  =  z-f  iy  =  0(t*) 
as  t  -*  0, 

\  WT  (  1  «n  ^/Ef^  ..  1 

+  ■  (35) 


2art 

d.  -  tdy  = 


This  result  can,  of  course,  be  obtained  directly  by  replacing  tanb(xz/4)  with  xz/4  in  (19);  the 
differential  equation  (24)  for  A  then  has  cosh  k  replaced  by  1  and  sinh  k  by  k. 

The  behavior  of  the  local  solution  obtained  as  r/f*  — » 0  is  more  complicated,  lb  expand  the 
integral  in  (35)  it  is  convenient  to  evaluate  the  following  two  integrals  separately; 


■'-Cl 

I,  ■ 


(«-<**  _i)+(i_!i^) 
in\^f 


k 


The  integral  /]  contains  a  term  that  cancels  the  logarithmic  term  in  (35).  In  Jj,  repeated  integra¬ 
tion  by  parts  provides  a  power  series  in  r/f  *.  Another  kind  of  term  iq>pear8  when  y/k  t  and  kx  are 
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both  large  aod  of  the  aame  order.  Thia  contributioo  ia  found  by  changing  the  integration  contour 
and  choosing  the  path  of  steepest  descent  from  k  =  t^/(Az^).  Substituting  these  approximations 
into  (35)  finally  gives 


for  t  — *  0  and  z/t’  — »  0,  where  7  =  0.577---  is  the  Euler  constant.  The  corresponding  local 
solution  for  the  surface  elevation  is 

1 = -^  [(i»  X + ■' - ’> + '(?) + 2(?)’ + •  ■•]  - -  f ) + •  •  •  +  o(«’)- 


This  result  is  identical  to  that  of  Roberts  (1987)  except  that  the  term  ln(irt^/4)  replaces  Int^. 
This  discrepancy  is  due  to  the  difference  in  the  problems:  Roberts  (1987)  considered  an  infinitely 
deep  fluid  with  a  finite-depth  wavemaker,  —  1  <  y  <  0 

One  important  implication  of  (32)  concerns  the  behavior  of  the  dynamic  contact  angle  between 
the  free  surface  and  the  moving  wall.  This  is  of  particular  interest  for  realistically  small  values 
of  the  nondimensional  surface  tension  T.  As  T  — »  0  it  is  obvious  from  (20)  that  a  singular- 
perturbation  problem  arises  near  the  contact  line.  The  sc^utira  to  this  problem  can  be  recovered 
from  (32).  Differentiation  of  (32)  with  res(>ect  to  x  and  introduction  of  the  transformation  k  =  kx 
yields  _ 

where  the  tanhk  factor  has  been  replaced  by  1  for  small  x,  as  in  (34).  This  form  is  useful  for 
small  T,  in  particular  for  T  -*  0  and  =  0{T^f*).  Now  if  */t*  -►  0  with  xfT^f^  fixed,  the 
cosine  term  in  the  mtegrand  oscillates  rapidly,  and  the  contribution  of  this  term  to  the  integral 
is  small.  The  remaining  term  can  be  integrated  to  give 


sin  kdk 
k[\+T'^)' 


»j,  =  (a  -  «c)  exp(-^)  ~a  +  0(a*), 


(38) 
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as  zfi^  — »  0  and  T  — »  0  with  zfT^I^  fixed.  Thus,  the  dynamic  contact  angle  approaches  the  static 
contact  angle  in  the  limit  of  zero  wavemaker  acceleration.  However,  when  the  surface  tension  is 
zero,  (38)  beccxnes  tjr  =  —a  as  z  — >  0  in  agreement  with  (37).  This  also  agrees  with  Roberts’ 
(1987)  findings  and  implies  a  jump  in  the  contact  angle  at  t  =  0  if  T  is  neglected.  For  small  but 
nonzero  T,  the  discrepancy  is  explained  by  the  large  curvature  —  Kjy/T  at  z  =  0,  as  shown 
by  (38). 

The  integral  in  (32)  is  evaluated  numerically  as  the  sum  of  integrals  from  zero  to  M  and 
M  to  infinity.  The  value  for  A/  is  chosen  to  be  as  large  as  10*  for  small  z  and  i,  and  as  small 
as  10~*  when  either  z  or  t  is  large.  The  integral  from  zero  to  Af  is  evaluated  using  10-point 
Gauss-Legendre  and  21-point  Kronrod  formulas  on  both  halves  of  the  adaptive  subintervals.  The 
selection  of  the  subinterval  is  based  on  the  maximum  absolute  error  estinute  of  10~^.  Due  to 
the  rapid  oscillation  cS  the  integrand,  the  integral  from  M  to  infinity  is  obtained  using  Filon’s 
method. 

Comparisons  of  (32)  with  the  small-tiine  s<4ution  and  the  local  solution  are  iUustrated  in 
Fig.  2.  For  small  time  (t  =  0.1)  and  T  s  0,  (32)  agrees  with  the  local  solution  (37)  near  the  wall 
(z  <  (’),  and  with  the  smaU-time  solution  (33)  sufficiently  far  frmn  the  wall.  Fig.  3  shows  the 
free-surface  configuration  at  small  time  for  two  different  ^ales  when  k  is  zero.  A  numerical  value 
of  the  nondimensional  surface  tension,  T,  few  pure  water  at  20*C  with  an  undisturbed  depth  of 
1  m  is  about  0.74  x  10~*,  represented  by  T  =  10~*  in  Fig.  3.  When  the  surface  tension  is  zero, 
small-scale  waves  (or  wiggles)  can  be  observed  very  near  the  wall,  as  also  noted  by  Roberts  (1987). 
These  wiggles  are  suppressed  in  the  presence  of  surface  tension,  in  which  case  the  static  contact 
angle  (90*  in  this  case)  is  retained.  At  a  given  time,  surface  tension  also  decreases  the  contact-line 
elevation.  Far  from  the  wavemaker  the  effect  of  surface  tension  becomes  less  important.  These 
effects  become  more  obvious  as  we  proceed  to  a  step  velocity. 
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4  Step  velocity 


The  wavemaker  velocity  given  by  (3)  when  9  s  0  is  a  step  velocity,  i.e.,  the  wavemaker  initially 
at  rest  is  suddenly  set  in  motion  at  t  =  O'**  with  a  constant  velocity. 

When  the  velocity  of  the  wall  C  is  small  compared  with  y/gK,  the  expansions  in  (QH^O) 
for  small  FVoude  number  can  be  used  as  before.  This  time,  the  Fteude  number  a  and  the 
dimensionless  velocity  u(f)  are  given  as 


Of  = 


and 


u{t)  =  1. 


With  u(t)  =  1,  the  solution  to  (24)  with  the  initial  conditions  in  (27)  is 


(39) 


(40) 


From  (16)  and  (22), 


=253  n  (I -CM  fit) 

*5  ^  Jo 


coshk(v+l)  .  dk 

— SSI— 


«a0  " 

After  cancellation  of  the  logarithmic  terms.  The  vertical  velocity  of  the  free  surface  becomes 

dk 


tanhkcoskx" 


(42) 


(43) 


- 1 

Equation  (19)  then  gives  the  free^urfaoe  elevation  for  a  =  0  as 

It  is  interestiiig  to  note  that,  sinoe  the  problem  is  linear,  (42)  and  (43)  can  be  derived  directly 
by  differentiating  (31)  and  (32)  with  respect  to  time.  This  procedure  is  equivalent  to  Roberts’ 
(1987)  use  of  a  convolution  integral. 

As  m  the  case  of  a  ramp  velocity,  (43)  with  7*  =  0  can  be  shown  to  agree  with 

t)  =  — to(tanh  ^)  +  0(at<*) 

sufficiently  far  from  the  wall,  for  O  t^,  and  with 

f/,  irt*  .  .  1  .  48crf.*../a  ,f’  x. 

+ •  •  ■ + 0(0’) 
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near  the  contact  line,  tot  x  <.t^.  These  equations,  the  small-tiine  sdution  and  the  local  solution, 
can  be  obtained  either  by  tedious  analysis  or  by  direct  time-differentiation  (33)  and  (37). 

Although  the  solution  procedure  differs  only  slightly  from  that  for  the  ramp  velocity,  subtleties 
arise  in  the  interpretation  of  (42).  If  we  hold  x  fixed  and  let  f  — »  0  in  (42),  and  then  let  z  — *  0, 
the  vertical  velocity  becomes  infinitely  large,  as  in  the  evaluation  of  (42)  for  z/f’  >  1.  If  instead 
we  reverse  the  order  in  which  these  two  limits  are  taken,  the  singular  behavior  disappears,  as  can 
be  seen  if  the  expansion  of  t;  for  z  C  is  substituted  in  (27).  The  nonuniformity  is  discussed 
below. 

If  (43)  is  rewritten  in  terms  of  k  =  kx,  it  can  be  shown  that  ty  =  Oiaf-^x)  as  t,z  — »  0 
with  x/t^  fixed.  This  behavior  of  q*  is  also  implied  by  the  expansions  for  large  and  small  z/t’. 
However,  since  the  assrmptotic  representaticm  requires  |i]b  |<^  I>  ^d  therefore  z  or^,  a  different 
inner  solution  is  needed  for  i  =  0{a),  x  —  0(a^).  The  conclusicm  can  also  be  anticipated 
from  dimensional  considerations.  The  relevant  parameters  for  points  near  the  contact  line  at 
small  time  should  be  C  and  g  rather  than  g  and  A,  so  that  the  proper  reference  length  and  time 
are  C^/g  =  at*h  and  C/g  =  at(fc/p)^/*,  respectively.  The  correspcmding  coordinates  are  t/a  and 
z/o*. 

The  correct  formulation  dt  an  inner  problem  for  small  t  and  z  can  be  shown  to  require  the 
full  nonlinear  free-surface  conditions.  The  proper  asjrmptotic  form  is  inferred  from  the  condition 
that  the  inner  solutions  for  q,  match  with  expansions  of  the  previous  solutions  obtained 

from  (41)  and  (43)  as  z,y,f  — *  0.  When  (43)  is  replaced  by  a  representation  analogous  to  (34), 
it  is  seen  that  the  largest  term  in  the  inner  solutkm  for  q  must  match  with  a  term  0{od  In  a). 
The  corresponding  y-coordinate  should  be  measured  from  this  first  approximation  to  the  surface 
elevation.  Matching  the  velocitieo  implies  that  —  O(o)  and  dy  =  0(a  In  a)  in  the  small-scale 
sdution. 

The  above  considerations  suggest  a  solution  ot  the  form 

d  =  (o*ln’a)di(x,y,<)-l-('?’*!ss)^(z,y,f)  +  o*ds(*,y,i)+  ■■  (44) 
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II  =  (a*lno)^(*,<)  +  o’%(x,f)+*- 


(45) 


where 


.  _  Z  y-(Q,»btt)^|(z.<)  .  _  £ 

*  ~  /yS  ’  *  “  ’  ft  ■ 


ft^  ■  ■  a*  ft 

Terme  O(o*  In^  o)  in  the  dynamic  free-surfsee  condition  (5)  give  +4if  =  0  and  no  =  di(0- 
Matching  gives  ^  — »  0  and  ^  — »  2/z  as  z  — »  oo;  since  the  complex  velocity  ^  is 

bounded  everywhere,  it  is  constant.  Terms  0(a*  In^  a)  in  (5)  then  give  4ii  —  ~ 

where  ~  2/*,  and  so  qj}  depends  only  on  t;  the  kinematic  condition  (4)  shows  that  =  2/t. 
Finally,  it  follows  from  (5)  that  —  ~b-  The  expansions  (44)  and  (45)  therefore  become 


^  =  (ft*  In’ at)'^  +  (of*toft)(~-^)  +  ft*^(z,y,  <)  +  ••• 

ij  =  (ft*lnft)“ +  ft’fe(z,  £)  +  ••• 

Conditions  to  be  satisfied  by  ^  and  %  at  y  =  ^  are  determined  from  (4)  and  (5): 

+  ^)  =  0 

^  ^  —  ^ih»  —  0. 


(46) 

(47) 

(48) 

(49) 


That  is,  the  full  nonlinear  freeeurface  conditions  are  required  for  ^  and  ^  and  are  to  be 
evaluated  at  the  unknown  locatkm  y  =  The  condition  (2)  at  the  wavemal»r  leads  to 


and  is  to  be  evaluated  at  the  actual  location  of  the  wall  z  =  f .  Thus,  the  linearised  formulation 
for  small  Fkoude  number  a  fails  when  t  s  O(o)  and  z  =  0(ft*).  Here,  a  full  notlinear  problem 
nnist  be  solved,  with  some  added  terms  invdving  In  at. 

The  large-time  behavim  of  the  fluid  motion  is  also  of  interest.  The  free-surface  configuration 
for  large  tinoe,  but  still  I  <  (l/o),  cno  be  obtained  by  an  asymptotic  evaluation  of  (43).  If  q  =  0 
initiaUy  (z  =  0),  the  free-surface  elevation  (43)  can  be  written  as 
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Ab  t  — »  oo,  the  largest  contribution  to  the  integrals  occurs  for  smaU  k,  at  k  s  0{l/(t  —  x)}.  If 
X  is  not  close  to  t,  it  is  sufficient  for  a  first  approximation  to  replace  ^  by  it  in  the  integrand 
and  to  omit  Tk^  in  the  denominator.  The  result  equals  a  for  x  <  f  and  zero  for  x  >  f  (recall 
that  dz/dt  =  1  corresponds  to  the  speed  y/pi  of  a  shaUow-water  wave).  This  approximation, 
however,  neglects  (k  —  0)t  ^  it  |t  —  x |  and  therefore  fails  near  the  wave  front,  z  =  t;  i.e.,  since 
i  —  /9  -x  (1  —  3T)k^/6  and  k  =  0{l/(i  —  x)},  it  has  been  assumed  that  |  f  —  x  (*<;  t.  When 
I  —  z  =  0(f^^’),  the  cubic  term  in  k  must  be  retained,  and  the  surface  elevation  becomes 


i,  =  ojf“Ai(Ode+  •.  (51) 

where  A  =  (1  -  3r)"‘/*(2/f  )‘^*(x  —  f)  and  Ai  is  the  Airy  function.  For  A  >  0  and  x  -  f  »  f 
17  is  exponentially  somU;  for  A  <  0  and  <  f  —  x  <  f, 

=  (‘'t)  “•  alrrjfj  '('•7)  +4  +  ••••  <“> 


Wh«u  f  —  X  =  0(<)  behind  the  wave  fr<»t,  the  first  integral  in  (50)  has  a  stationary  point  at 
k  =  0(1)  and  contributes  a  term  0(at~^f^)  that  matches  with  (52).  In  the  secmd  integral  there 
is  no  stationary  point  when  x  >  0,  and  the  integration  contour  can  be  deformed  to  lie  somewhat 
away  from  the  real  axis  in  the  complex  k>plane.  The  largest  contribution  is  near  k  =  0,  giving 
the  value  x/2  with  exponentially  small  error. 

Therefore,  as  time  becomes  large,  the  contact-line  elevation  ^proaches  a  value  equal  to  the 
Fkoude  nundber,  and  behind  the  srave  front  the  free  surface  can  be  approximated  by  a  wavetrain 
superimposed  on  a  flat  surface  of  the  same  height  as  the  contact  line.  The  amplitude  of  this 
wavetrain  is  0(<rt~'/*),  increasing  to  0(a)  near  the  wave  front.  Surf .  .  tension  increases  the 
frequency  and  decreases  the  amplitude  decay  rate  of  the  wavetrain.  Beyond  the  wave  front  the 
free-surface  elevation  decreases  exponentially  to  the  undisturbed  value  of  zero.  The  width  of  the 
wave  front  increases  in  proporticm  to  so  that  the  slope  ot  the  free  surface  near  x  =  t  is 
0(ot-‘/»). 


In  Fig.  4,  the  free-surface  elevation  is  shown  at  several  different  times  for  zero  initial  elevation 
(k  =  0)  and  Cor  three  values  of  the  nondimensional  surface  tension.  Near  the  wavemaket  (x  <;  f*), 


the  free  surface  m  made  up  of  ao  infinite  number  of  wiggles,  which  can  be  approximated  by  a 
local  solution.  Surface  tension  suppiuaca  these  wiggles,  but  the  effects  of  surface  tension  are  seen 
to  decrease  as  the  distance  from  the  wavemaker  increases.  This  explains  why  the  agreement  with 
experimental  measurements  (Dommermuth  et  al.,  1988)  is  good  even  though  surface  tension  is 
neglected  at  moderate  distances  from  the  wavemaker.  Also  shown  in  Fig.  4e  is  the  asymptotic 
solution  (51)  evaluated  using  the  Taylor  series  and  the  asymptotic  representations  in  Abramowits 
ic  Stegun  (1985),  which  is  in  good  agreement  with  the  numerical  evaluation  of  the  exact  solution 
(43)  near  the  wave  front. 

The  distance  from  the  wavemaker  occupied  by  the  wiggles  increases  with  time  in  agreement 
with  the  local  solution.  The  contact>iine  elevatkm  increases  to  a  maximum,  and  then  oscillates 
to  converge  to  the  fVoude  number  at  large  time  (Fig.  5).  Fig.  4d  shows  that  for  certain  values 
of  time  the  surface  tension  actually  makes  the  contact-line  elevation  higher,  as  also  shown  in 
Fig.  5.  The  amplitude  and  frequency  of  the  wiggles  near  the  wavemaker  decrease  with  time 
until  the  free  surface  becomes  flat,  as  indicated  by  the  analysis  (50V(52).  The  extent  of  this  flat 
region  increases  with  time,  so  that  in  the  limit  mi  —*  oo  the  behavior  of  the  free  surface  can  be 
^proximated  by  a  hydraulic  jump,  with  the  jump  location  at  x  =  f .  This  corresponds  to  the 
phase  velocity  in  shallow  water,  as  can  also  be  semi  by  iropoamg  conservation  laws  in  a  simple 
control-volume  analysis. 

5  Exponential  wavemaker  velocity 

As  a  general  example  that  includes  the  step  and  ramp  velocities  as  limiting  cases,  we  consider 
a  wavemaker  velocity  that  has  a  finite  jump  in  acceleration  at  f  =  0  and  approaches  a  constant 
value  I/q  as  t  — » oo.  The  exponential  form 

U(t)  =  Uo 

where  r  is  a  characteristie  time,  exhibits  this  behavior.  The  limits  r  —>  0  and  r  -» oo  correspond 
to  the  step  velocity  and  the  ramp  velocity,  respectively.  The  expansions  for  small  FVoude  number, 
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(9)-(10),  are  applied  with 


and  t»(t)  =  1  —  e~**,  (54) 

as  the  Froude  nund>er  and  nondimensional  velocity,  respectively.  Here  b  —  \/h/(gT^).  Now,  (24) 
and  (27)  give 


-  "7 (^1  -  cos^  -  0  j  . 


Equations  (16)  and  (22)  then  give  the  first-order  velocity  potential  as 


di  =  2(l-e-“)5;^e 
n«0*“ 


1 


Sltt 


—CO»0t  — 


0 


0t~**  +bmi0t-  0C(»0t\  dk 

l^  +  0»  ) 


(55) 


(56) 


The  vertical  vdocity  on  the  free  surface,  after  cancellation  of  the  singular  terms,  beccmes 


‘■■ir 


cos^tanhhcoshx 


r 


il*t~**  -  0bBm0t  + 0^  cc»0t  ^  ^  dk 

- - tanhtcost,-, 


(57) 


and  the  free-surface  elevation  becomes 


_2a 

~  ^  Jo 


0m0t 


k^l  +  Tk*) 


cmkxdk+ 


2a  f*"be~*'-bem0t-0m0t^  ^  dk  , 

-j,  - 


(58) 


With  t  fixed,  if  6  — »  0,  (56)  approaches  the  sdution  (32)  for  the  ramp  velocity  with  a  replaced 
by  a6;  if  5  -*  oo,  it  approadies  the  stdutkm  (43)  for  the  step  velocity.  The  first  term  in  (58)  is 
identical  to  that  for  the  step  velocity,  and  the  other  term  decays  to  sero  as  <  — » oo.  Consequently, 
as  t  oo  the  bdiavior  of  the  fluid  eventually  follows  that  of  the  st^  velocity  regardless  of  the 
startup  process.  As  t  -*  0  a  qf>ecial  ease  arises  when  b  -*  oo  such  that  bt  becomes  constant  in 


the  limit — the  solutkai  then  depends  oo  the  value  ofM.  IfM— >0otM-»oo,  asf— >0,  the 
mall-time  solution  is  recovered  for  the  ramp  and  the  step  velocity,  respectively. 

We  have  shown  in  the  previous  sections  that  the  contact  an^  of  an  invisdd  fluid  remains 
unchanged  from  its  initial  statk  state  in  the  presence  of  surface  ieuskm.  We  now  examine  the 
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relationahip  between  the  Kceler«ti<»  of  the  wnvemnker  and  the  djmamic  contact  angle  when  the 
capillary  effect  is  absent. 

When  the  surface  tension  is  sero  (T  =  0),  differentiation  of  (58)  with  respect  to  x  and  the 
transformation  k  =  kx  yield 


»?.  =  - 


—  r  (l - \  Jitanh-sin(N/itanhi-~)8iniS 

»  io  V  ei’  +  ttanhi/Y  *  >/*  * 

[e“*‘  -  cos(t/fctanhT-^)] - * — rsinhdh  +  0(a^). 

*  JO  I  J  4-  k  tanh  ^ 


(59) 


When  «  -» 0  and  x/i^  — » 0,  this  can  be  evaluated  as 


ij*  =  — a6e~**  H - 


(60) 


Therefore,  the  slope  of  the  free  surface  very  near  the  contact  line  jumps  instantaneously  to  a 
finite  value,  which  depends  on  the  Ftoude  numb^  and  the  exponent  h,  and  decays  to  sero  again 
as  tune  becomes  large.  Since  the  decay  rate  increases  with  6,  it  is  not  surprising  to  observe  tie 
90*  contact  angle  for  the  step  velodty  even  at  an  inflnitesima]  time  (Fig.  4).  Since  the  limit  5  — » 0 
corresponds  to  the  ramp  velocity  for  finite  time,  (60)  can  easily  be  shown  to  be  consistent  with 
the  previous  result,  if  the  difference  in  definitioo  of  the  Fteude  numbers  is  taken  into  account. 

As  shosm  in  the  previous  section,  the  linear  scdution  for  the  step  velocity  is  not  valid  for  very 
small  time,  and  a  fully  nonlinear  formulation  is  required.  This  can  be  more  easily  understood  by 
examining  the  limitation  to  be  imposed  on  the  time  constant  (or  b)  for  the  present  wavemaker 
if  the  linear  solution  is  to  be  unifonnly  valid.  As  we  proceed  to  the  next  order,  O(a’),  terms 
proportional  to  ei*be~**  wiU  appear.  Then,  for  5  >  1  the  expansioos  (9)-(10)  are  valid  fw  all 
time,  including  t  (1/5)  only  when  b  <  (1/a).  If  5  =  0{l/a),  nonlinear  terms  are  required 
when  bt  s  0(1).  Therefore,  the  nonlinear  effects  cannot  be  neglected  for  a  nqiidly  accelerated 
wavemaker  (r  =  0{V,/t)).  This  is  also  cmsistent  with  (60)  when  it  is  combined  with  the 
kmematk  boundary  condition  on  the  free  sur&ce  (4),  which  gives  the  same  criteria  for  the 
validity  of  linearisation  as  above. 

The  free^urfaee  configurations  for  varioas  vahiee  of  5  are  shown  in  Fig.  6  whoi  surfoce  toi* 
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sion  and  initial  free-surface  elevation  are  absent.  They  resemble  those  for  the  ramp  velocity 
when  b  is  small,  and  those  for  the  step  velocity  when  b  is  large.  As  the  acceleration  of  the  wave- 
maker  increases  (as  b  becomes  larger),  the  amplitude  of  the  wiggles  grows  and  the  contact  angle 
approaches  90”,  which  is  consistent  with  the  above  analysis. 

6  Harmonic  wavemaker  velocity 


Simple-harmonic  motion 


A  wavemaker  motion  of  greater  practical  interest  is  the  periodic  oscillation,  for  which  the 
velocity  of  the  wall  is  given  as 


U(t)  s  f/baio 


t  >  0, 


(61) 


where  Uo  is  the  maximum  velocity  of  the  wall  and  0  is  the  frequency  of  the  oscillation.  The 
well-known  linear  steady-state  solution  to  this  problem  was  obtained  by  Havelock  (1929)  and  has 
been  extended  by  many  others.  The  Fourier-integral  method  adopted  in  this  study  will  lead  to 
a  transient  solution  that  agrees  with  the  steady-state  solution  as  f  — » oo. 

The  FYoude  number  is  defined  as  in  the  previous  section,  and  the  normalisation  of  the  wave- 
maker  velocity  gives 


«(()  =  sinwf. 


(62) 


where  the  nondimensional  frequency  of  the  wavemaker  oscillation  is  w  =  Qy/h/g,  so  the  solution 
for  (24)  that  satisfies  (27)  is  now 


^  0  Pmaut  —wnn0t 


(63) 


Therefore,  the  complete  first-order  velocity  potential  is 

=2smwf  V-i-e-*-»sink,.yf  1  T r  k(y 1)  dk 

(u»a-/?*)coshfc  cos»*^,.  (64; 

Again,  the  singular  terms  in  the  vertical  velocity  are  cancelled  to  give 

j  2w  umnut- 0m0t  ^  ,  dk 

- tanh*  costs—  (65) 
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on  the  free  turfaee.  The  free-nirfaee  elevation  is  then 


2au  coewt  —  coe^  ,  „  3, 


Far  from  the  wavemaker,  the  asymptotic  evaluation  of  (66)  is  possible  in  the  limit  as  t  — »  00. 
As  X  and  t  become  large,  the  largest  contribution  to  the  integral  in  (66)  occurs  in  the  neighborhood 
of  k  =  ko,  where  ho  is  the  positive  real  root  of 

ko(l  +  rkg)  tanh  ho  =  (67) 

Since  this  is  just  the  dispersion  relation,  ko  is  the  wavenumber  that  would  be  observed  for  waves 
with  the  single  frequency  w.  Therefore,  (66)  can  be  written  as 


2au 

*■  Ju-* 


’•**  coBut  cox kx  -  \  cos(t«  -f  ^)  —  ■^  cos(hx  —  0t)  dk 


where  e  is  a  smaD  number  such  that  1  <  (l/<)  <  x.  The  integrand  in  (68)  can  be  expanded 
about  k  =  ko  and  the  resulting  equation  can  be  easily  simplified  after  setting  {k  -  ko)*  =  k: 


ataohko 


-ff- 

i  Jo 


2coowt  mnkoX’ 


to..) -.^0  Ji-f  -  ■  ■ .  (69) 

Here,  Cg  is  the  group  velocity  of  the  gravity-c^illary  wave  with  wave  number  ko  and  is  given  by 


When  X  >  Cft,  the  approodmation  (69)  is  sero;  when  x  <  Cgt,  (69)  gives 

o tanh ko  .  ..  . 

»»* - j^j^^sin(to*- «<)  +  •••.  (71) 

which  describes  the  free-surface  configuration  in  a  region  behind  the  wave  front  but  far  ahead  of 
the  wavemaker.  In  the  absence  of  surface  tension  (T  s  0),  the  approodmation  (71)  agrees  with 
the  steady*state  solution  sway  from  the  wavemaker  obtained  by  Havelock  (1929).  As  for  the  step 
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velocity,  the  behavior  of  the  free  surface  near  the  wave  front  (z  =  Cgt)  could  be  obtained  by 
extending  the  above  analysis  to  higher  orders. 

The  behavior  of  the  contact  angle  in  the  absence  of  surface  tension  is  obtained  by  differenti¬ 
ating  (66)  with  respect  to  x  and  transforming  h  to  £  as  before,  which  gives 


cosw<-cos(^£tanh^^] 


This  can  be  simplified  to 


Jig  s  —au  ct»ut  +  •  •  • 


(72) 


as  x/t*  — »  0  and  z  — •  0.  The  discontinuity  in  the  contact  angle  at  t  =  0  is  consistent  with  the 
previous  results.  It  can  be  easily  observed  from  (72)  that  the  contact  angle  oscillates  with  a  90* 
phase  shift  from  the  wavemaker  velocity. 

Fig.  7  shows  the  fiee-surface  configuration  at  large  time  for  two  different  wavemaker  frequen¬ 
cies.  The  harmonic  wavetrain  of  Havelock  (1920)  is  observed  between  the  wavemaker  and  the 
harmonk-wave  front,  z  s  Cgi.  The  amplitude  and  frequency  of  this  wavetrain  can  be  obtained 
easily  from  (67)  and  (70)-(71).  As  the  undisturbed  free  surface  is  approached,  ahead  of  the 
harmook-wave  front,  we  obanve  a  second  wave  front  that  travels  at  the  maxunum  phase  velocity 
dx/dt  sz  1,  the  value  for  shallow-water  waves.  The  waves  between  these  wave  fronts  are  seen 
to  have  decreasing  amplitude  and  wavenumber.  The  largest  change  of  amplitude  occurs  near 
z  =  Cgt  in  both  Fig.  7a  and  7b. 


Harmonic  analysu  of  a  general  wavemaker  velocity  ’ 


As  a  final  example,  the  velocity  of  a  wavemaker  given  by  a  Fourier  cosine  series 

/f 


«(0  =  52 «n  COS(w,t  -  ^n) 


(73) 


nml 


u  considered.  A  straightf<»ward  applicatkm  of  (27)  gives 


A  =  i 


vPmSlt  ^  fw„ sinfn sin^  -  fiet»0„«oe^  +  0cax(unt  - f«)] .  (74) 

nal  ^ 
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We  then  follow  the  tame  procedure  ac  before  to  obtain 

Wfi linfwwt  —  dn)  -f  fa>n aiP 9nCoa$t  —  0e(m6n*oi0t 


=  2o 

'  Jo  [^l 


w; 


dk 


tanhh  coat*— +O(o*). 

(75) 


We  take  N  ^72  and  T  =  0  and  uae  the  Fourier  come  coefficienta  of  the  wavemaker  velocity 
provided  by  Dommermuth  et  al.  (1988).  Fig.  8  coii4>area  the  free-aurface  elevationa  againat 


time  at  three  different  locationa  to  the  reaulta  Dommermuth  et  al.  (1988)  obtained  by  the 
boundary-ini^ral  method  with  lineariied  free-autface  conditiona.  We  have  choaen  to  include  only 
the  extrema  of  their  figurea  for  clarity.  Alao  ahown,  when  aufficiently  different  from  the  linear 
computational  are  their  wave-probe  meaaurementa.  For  a  moderate  diatance  from  the  wavemaker 
(x=3.17),  the  free-aurface  elevation  (75)  ahowa  good  agreement  with  the  experimental  meaaure¬ 
menta.  Farther  away  from  the  wavemaker  (x=9.17),  the  nonlinear  effecta  have  accumulated,  and 
the  agreement  becomea  leaa  aatkfiactory.  In  aD  caaea,  very  cloae  agreement  in  obaerved  between 
the  preaent  Fourier-integral  aolution  and  the  numerical  aolution  ot  Dommermuth  et  al.  (1988). 
Since  the  linear  aolutko  (75)  ia  exact  and  can  be  evaluated  aa  accurately  aa  we  chooae,  the 
amall  difference  between  (75)  and  the  numerical  aolution  can  be  attributed  to  the  difficultiea  in 
the  boundary-integral  computation  anaodated  with  the  contact  line.  Aa  ahown  in  the  previoua 
aectiona,  aurfaee  tenaion  aflecta  the  free-aurface  eonfiguratioa  primarily  very  near  the  wavemaker 
(*  <  f’)  and  tot  amall  time,  ao  it  in  neglected  in  tbeae  compariaona. 


7  Concluding  remarks 


Ib  avoid  aa  artificial  aingularity  at  the  contnet  line  between  the  free  aurfaee  and  the  waven»ker 
introduced  by  the  mull-time  expaaaioo,  a  Fourier^integral  method  ia  devek^wd  for  amall  Fteude 
nunober.  Thia  method  yielda  aolutioaa  that  are  uaifonnly  valid  for  general  wavemaker  velocitiea 
that  need  not  be  given  aa  poweta  of  time.  It  alao  allowa  the  atndy  of  the  capillary  effecta  with 
little  added  effort. 

In  the  abaeace  of  anrfoee  tenaion,  an  infinite  number  of  amaO-acale  wig^  ie  preaent  near  the 
wavemaker,  aa  ahown  alao  in  the  local  aohitioo  of  Roberta  (1967)  for  amall  time,  and  the  cmtact 
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angle  hai  a  junq>  at  t  s  0.  Surface  tenaion  suppresaea  the  wigglea  and  maintaina  the  contact 
angle  at  ita  initial  atatic  value.  Far  from  the  wavemaker,  eflecta  of  aurface  tenaion  become  leaa 
important.  For  connatent  and  realiatic  reaulta,  aurface  tenaion  ahould  be  considered  near  the 
contact  line. 

When  the  acceleration  of  the  wavemaker  is  aufficiently  large,  the  present  linear  solution  is 
not  uniformly  valid  near  the  wavemaker  few  very  small  time.  A  correct  inner  solution  for  these 
conditions  requires  a  full  nonlinear  formulation. 

The  large-time  behavior  for  the  wavemaker  moving  with  constant  velocity  is  also  obtained. 
The  contact-line  elevation  approaches  a  value  equal  to  the  FVoude  number,  and  the  free  surface 
behind  the  wave  front  can  be  approximated  by  a  wavetrain  superimposed  on  a  flat  surface.  Be¬ 
yond  the  wave  front,  which  nooves  with  the  phase  velocity  few  shallow  water,  the  free-surface 
elevation  decreases  exponentially  to  the  undisturbed  value,  aero.  For  the  simple-harmonic  wave- 
maker,  the  large-time  behavior  agrees  srith  the  steady-state  solution  of  Havelock  (1929)  behind 
the  wave  front  but  frur  from  the  wavemaker. 

A  general  wavemaker  velocity  given  by  a  Fourier  cosine  scries  is  considered,  and  the  free- 
surface  elevation  is  compared  with  the  computatiooal  and  experimental  results  of  Dommermuth 
et  al.  (1988).  The  present  Fourier-integral  solution  agrees  very  well  with  the  linear  numerical 
results.  The  agreement  with  the  wave-probe  measurements  is  excellent  at  moderate  distances 
from  the  wavemaker  for  all  time  and  becomes  less  satisfactory  farther  aw^r  from  the  wavemaker 
for  large  time,  when  the  nonlinear  effects  have  accumulated. 
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LIST  OF  FIGURES 


Figure  1.  Wavemaker  configuration. 

Figure  2.  Ftee-surface  elevations  for  the  ramp  velocity  at  t  =  0.1  according  to  present  method 
( - )  and  compared  to  the  small-time  solution  ( - )  and  the  local  solution  ( . ). 

Figure  3.  FVee-surface  elevations  f<»  the  ramp  velocity  at  t  =  0.1  when  surface  tensicui  T  =  0 
( - ),  T  =  10-“  (-  -  -  -),  and  T  =  10"“  ( . ):  a)  large  x/t»;  b)  smaU  i/t*. 

Figure  4.  F^ee-surface  elevations  for  the  step  velocity  when  surface  tension  T 

T  =  10-“  ( - ),  and  r  =  10-“  ( . ):  a)  t  =  0.1;  b)  t  =  0.1,  x  <  t*;  c)  t 

X  C  t;  e)  t  =  100  with  the  asymptotic  solution  near  the  wave  front  ( —  •  — ). 

Figure  5.  Contact-line  elevation  for  the  step  velocity  when  surface  tension  T  =  0  ( - )  and 

T  =  10-“  ( - ). 

Figure  6.  FYee-surface  elevations  for  the  exponential  velocity  at  t  =  0.1  when  surface  tension 
T  =  0;  a)  6  =  1;  b)  6  =  10;  c)  6  =  50. 

Figure  7.  FYee-surface  elevations  for  the  simple-harnwnic  velocity  at  t  =  60  when  surface  tension 
T  =  0;  a)  w  =  1;  b)  w  =  2. 


=  0  ( - ), 

=  1;  d)  t  =  1, 


Figure  8.  FVee-surface  elevations  according  to  present  method  ( - )  and  compared  to  ex¬ 

perimental  measurements  (o)  and  computational  results  (•)  of  Dommermuth  et  al.  (1988):  a) 
X  =  3.17;  b)  X  =  5;  c)  x  =  10. 
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